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Abstract 

We construct a super symmetric (SUSY) SU(5) model with the flavor symmetry 
S4 x Z3 x Z4. Three generations of adjoint matter fields are introduced to generate 
the neutrino masses via the combined type I and type III see-saw mechanism. The 
first two generations of the the 10 dimensional representation in SU(5) are assigned 
to be a doublet of £4, the second family 10 is chose as the first component of the 
doublet, and the first family as the second component. Tri-bimaximal mixing in the 
neutrino sector is predicted exactly at leading order, charged lepton mixing leads 
to small deviation from the tri-bimaximal mixing pattern. Subleading contributions 
introduce corrections of order to all three lepton mixing angles. The model also 
reproduces a realistic pattern of quark and charged lepton masses and quark mixings. 
The phenomenological implications of the model are analyzed in datail. 
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1 Introduction 



So far there is convincing evidence that the so-called solar and atmospheric anomaly 
can be well explained through the neutrino oscillation. The mass square differences Am^,, 
Am^ tm and the mixing angles have been measured with good accuracy [TH3]. Global fit to 
the current neutrino oscillation data demonstrates that the observed lepton mixing matrix 
is remarkably compatible with the tri-bimaximal (TB) mixing pattern [I], which suggests 
the following values of the mixing angles 

sin 2 6> 12 = -, sin 2 6» 23 = -, sin0 13 = O (1) 

This simple structure of the mixing matrix suggests that there may be some symmetry 
underlying the lepton sector. Recent study [5T-TT9] (also in the context of grand unified 
theories [20j[2T]) showed that the discrete group A 4 is particularly suitable to produce the 
TB mixing at leading order (LO), if it is properly managed to be broken differently in the 
neutrino and the charged lepton sector. In the A4 based models, it seems very difficult and 
unnatural to generate the correct mass hierarchies and mixing for quarks. An interesting 
solution is to enlarge the symmetry group A4, two non-Abelian discrete group T' [22] and 
64 [23H30] have been investigated. We note that both groups have a doublet representation, 
which can be utilized to give the 2+1 representation assignments for the quarks. In the 
context of U(2) flavor group, this assignment has been known to give realistic quark 
mixing matrix and mass hierarchy [31]. The irreducible representations of T' are those of 
A4 plus three two dimensional representations 2, 2' and 2" with the multiplication rules 
2 © 2 = 2' © 2" = 3 © 1, 2 © 2' = 2" © 2" = 3 © V and 2 <g> 2" = 2' <g> 2' = 3 © 1", 
these ingredients allow us to reproduce the successful U(2) predictions in the quark sector. 
By working only with the triplet and singlet representations, T' is indistinguishable from 
A4, thus we can replicate with T' the successful construction realized within A± in the 
lepton sector. S4 is claimed to be the minimal group which can predict the TB mixing 
in a natural way, namely without ad hoc assumptions, from the group theory point of 
view [32] . Actually the exact TB mixing can be realized in the S4 flavor model [2"3ll2"l] . 
Moreover, the group S4 as a flavor symmetry, as is shown for example in Ref . [27H30] , can 
also give a successful description of the quark and lepton masses and mixing angles within 
the framework of SU(5) or 5*0(10) grand unified theory (GUT). For a review of discrete 
flavor symmetry models, please see the Refs. [3^135] . 

The SU(5) GUT is the simplest grand unified theory [33], in this case each generation 
of the standard model fields resides in 5 and 10 dimensional representations. To be spe- 
cific, one family of right-handed down quarks and left-handed leptons are unified in a 5 
and the rest fields of the family are in a 10. It is well-known that neutrino masses are zero 
at renormalizable level in the minimal SU(5). In GUT neutrino masses come naturally 
through the see-saw mechanism, where integrating out large masses leads to the appear- 
ance of small masses. However, this requires some extra matter fields or Higgs to be added 
below the GUT scale. A popular choice is to add at least two right-handed neutrinos which 
are SU (5) singlets, neutrino masses are generated through type I see-saw mechanism. The 
second choice is to introduce a 15-plet of Higgs, this is the SU(5) implementation of the 
so-called type II see-saw mechanism. The third choice is to generate neutrino masses 
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through the combination of type I see-saw with type III see-saw mechanism by adding at 
least one matter multiplet in the adjoint 24 representation [36H38] . 

In this work, we shall construct a supersymmetric SU(5) model, and the flavor symme- 
try group is S4 x Z3 x Z4, where the auxiliary symmetry Z3 x Z4 plays an important role in 
eliminating unwanted couplings, ensuring the needed vacuum alignment and reproducing 
the observed charged fermion mass hierarchies. We will introduce three generations of 
adjoint matter fields to generate the neutrino masses. We remark that some variants of 
SU(5) x S4 flavor models have been proposed in Refs. [28l[29], where three right handed 
neutrinos are introduced and the neutrino masses are generated via the type I see-saw 
mechanism. 

The paper is organized as follows. In section 2 we discuss the structure of the model, 
and the leading order (LO) results for fermion masses and mixings are presented. In section 
3 we justify the choice of the vacuum configuration assumed in the previous section, by 
minimizing the scalar potential of the theory in the supersymmetric limit. In section 4 
the subleading corrections to the vacuum alignment and the LO results of fermion masses 
and flavor mixings are discussed. In section 5 we study the phenomenological predictions 
of the model in detail. Finally, section 6 is devoted to our conclusion. 

2 The structure of the model 

In the following we present our 577(5) model in the framework of supersymmetry, which 
simplifies the minimization of the scalar potential greatly. The S4 group acts as a flavor 
symmetry of our model, the group S4 has already been studied in literature [311110] , but 
with different aims and different results. 5*4 is the permutation group of four objects, it has 
five irreducible representations li, I2, 2, 3i and 32, the group theory of 64 is presented in 
Appendix A. In addition to 5 matter fields denoted by F and the tenplet 10 dimensional 
matter fields denoted by T 12 3 , we introduce the chiral superfields A in the adjoint 24 
representation [37]. In the Higgs sector we introduce if 24, H5 and H§ in order to break 
the gauge symmetry SU (5) into the standard model symmetry and subsequently into the 
residual S77(3) c x U{l) em . Moreover, if 45 and are introduced to avoid the wrong 
predictions Mj = Mg, where Md and Mi represent the mass matrix of down type quark 
and charged lepton respectively. As usual, the flavon fields are introduced to spontaneously 
break the S4 flavor symmetry properly. The transformation rules of the matter fields, Higgs 
fields and the flavon fields under S77(5), S4, Z 3 and Z4 are summarized in Tabled] The 
first and the second generation of the 10 dimensional representations are assigned to be 
a doublet of S4, and the third generation of 10 to li of S4. This assignment is indicated 
by the heaviness of the top quark. There is the freedom of choosing the first family or the 
second family as the first component of the S4 doublet. In this work, the second family 
10 is taken to be the first component of the doublet, and the first family as the second 
component. If we assign the first family 10 as the first component of the doublet and 
the second family as the second component as usual, the down quark and strange quark 
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masses would be of the same ordeio without fine tuning unless some special mechanisms 
are introduced such as Ref. [28]. Three generation of 5 fields F are assigned to 3i of S±, 
and three generations of adjoint matter fields A are also assigned to be 3i of S4. Fermion 
masses and mixings arise from the spontaneous breaking of the flavor symmetry by means 
of the flavon fields. In the following, we shall discuss the LO predictions for fermion masses 
and flavor mixings. For the time being we assume that the scalar components of the flavon 
fields acquire vacuum expectation values (VEV) according to the following scheme 




We will prove it to be a natural solution of the minimization of the scalar potential in 
section 3. Furthermore we take the VEVs (scaled by the cutoff A) v x /A, v v /A, v^/A, v v /A, 
va/A and v^/A to be of the same order of magnitude about 0(\ 2 C ) with A c ~ 0.22 being 
the Cabibbo angle, and we will parameterize the ratio VEV/A by the parameter e. This 
order of magnitude is indicated by the observed ratios of up quarks and down quarks and 
charged lepton masses, by the scale of the light neutrino masses and is also compatible 
with the current bounds on the deviations from TB mixing for the leptons. We note that 
the assumed size of the VEVs can be partially explained by the minimization of the scalar 
potential, as it will be clearer in the following. 

2.1 Neutrino sector 

The LO superpotential which contributes to the neutrino masses is given by 

w u = y v {FA) u Hs + Xi{AA) 3lX + X 2 (AA) 2 ip (3) 

1 li we take T\ and T 2 as the first and the second components of the S4 doublet, i.e., (Ti, T 2 ) T ~ 2, then 
(TF) 3l ~ (TiFa + T 2 F 3 ,T 1 F 3 + T 2 F 1 ,T 1 F 1 + T 2 F 2 ) T and (TF) 32 ~ {T\F 2 - TiF^TxF* ~ T 2 F ll T 1 F 1 - 
T 2 F 2 ) T . We see that T\F\ and T 2 F 2 , which are related to the down and strange quark masses respectively, 
appear simultaneously as the third component of both the combinations (TF) 3l and (TF)a 2 . The opera- 
tors TFHir and TFH-^ combining with the flavon fields or the composition of flavons, which transform as 
3i or 32, contribute to the first two families down quark and charged lepton masses after the £4 and GUT 
symmetry breaking. As a result, the down and strange quark messes would be of the same order except 
for the case that the (TF) 3l and (TF) 32 relevant contributions to down quark or strange quark cancel 
with each other. We notice that the same view has been put forward by Altarelli et al. [31]. Whereas 
if we choose {T^f - 2 as we proposed, then (TF) 3l - (T 2 F 2 + T^^Fs + T 1 F 1 ,T 2 F 1 + T^)? 
and (TF) 32 - (T 2 F 2 - TxF 3 , T 2 F 3 - TyFi^Fx - TiF 2 ) T , the degeneracy between first two families down 
quark (charged lepton) masses is dissolved. 
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Table 1: Fields and their transformation properties under the symmetry groups SU(5), S4, Z 3 and Z4, 
where w = e l2lT / 3 = (-1 + iy/S)/2. 

The adjoint matter field A decomposes under the standard model as A = (p 8 , p 3 , P(3,2), P(3,2) ; Po) 
(8, 1, 0)©(1, 3, 0)©(3, 2, -5/6)0(3, 2, 5/6)0(1, 1, 0), obviously there are both 377(2) triplet 
P3 and singlet po with hypercharge Y" = in the model. The neutrino masses are generated 
through the type I (mediated by the SU(2) singlet po of A) and type III (mediated by the 
SU(2) triplet p^ of A) see-saw mechanism. The Dirac mass matrices is obtained from the 
first term in Eq. (13|) , 
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where M% and are the Dirac mass matrices associated with p 3 and p respectively. 
The last two terms in Eq.(|3]) lead to the Majorana mass matrices of p% and po 
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It is notable that the Majorana mass matrices of P3 and po are exactly the same. As a 
result, the mass spectrums of p 3 and po are degenerate. This degeneracy is violated at 
NLO by the Higgs H 2i . We note the VEVs (x) and (<p) are invariant under the action of 
£4 elements TSTS 2 , TST and S 2 , consequently the £4 flavor symmetry is broken down 
to the Klein four subgroup in the neutrino sector. The light neutrino mass matrix is the 
sum of type I and type III see-saw contributions 
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where 

a = Xiv x , b = X 2 v v (7) 

The above light neutrino mass matrix M v is 2 -h- 3 invariant and it satisfies the magic 
symmetry (M„)n + (M„)i3 = {M v ) 22 + (M u ) 32 . Therefore it is exactly diagonalized by the 
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TB mixing matrix 



UjM u U v = diag(mi, m 2 ,m 3 ) 



where m li2j3 are the light neutrino masses, in unit of ^y„v£ they are 

1 

mi 



m 2 



m 3 



\3a ■ 
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2\b\ 



|3a + 6| 



(9) 



The neutrino mass spectrum can be normal hierarchy (NH) or inverted hierarchy (IH). 
The unitary matrix U v is given by 



U u = U TB diag(e' 
Utb is the well-known TB mixing matrix 
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The phases a±, a 2 and a 3 are 
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According to Eq.([9]), the light neutrino mass spectrum is directly related to the heavy 
neutrino (p 3 or po) masses, it is determined by only two complex parameters a and b. In 
the following, we shall follow the method of Ref. [17] to analyze the light neutrino mass 
spectrum in detail. For the sake of convenience, we define 
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Note that the parameter r is real and positive and the phase 9 is between and 2ir. We 
can express r and the phase 9 in term of light neutrino masses as follows 
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Experimentally, only the mass square differences have been measured. For normal (in- 
verted) hierarchy they are 
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Figure 1: cos 8 as a function of the lightest neutrino mass m; for both normal hierarchy and inverted 
hierarchy spectrum, where we take Am 2 o( = 7.67 x 10 -5 eV 2 and Am 2 tTO = 2.46(2.45) x 10 _3 eV 2 for 
normal hierarchy (inverted hierarchy). 



By expressing the neutrino masses in terms of the lightest neutrino mass mi (mi = mi for 
NH and m 3 for IH), Am^ and Am^ tm , cos 9 becomes a function of the lightest neutrino 
mass mi. We display cos# versus mi in Fig |T] for both normal hierarchy and inverted hi- 
erarchy spectrum. Imposing the condition | cos^l < 1, we obtain the following constraints 
on the lightest neutrino mass 

mi > O.OlleV, for normal hierarchy 

m 3 > 0.028eV, for inverted hierarchy (17) 

The results presented so far are of course approximate since the model gets corrections 
when higher dimensional operators are included in the Lagrangian. If the mixing of the 
left handed charged lepton is neglected, we can straightforwardly derive the neutrinoless 
double decay parameters \m ee \ 




mee | = _iW2--f + 2-f (18) 



2 ' 2 

m\ m\ 



2.2 Up quarks sector 

The LO superpotential invariant under symmetry group SU (5) x S 4 x Z 3 x Z 4 , which 
gives rise to the masses of the up type quarks after S4 and SU(5) symmetry breaking, is 
given by 



y t T 3 T 3 H 5 + %sTTO®H s + ^TT 3 (^) 2 if 5 + y -^TT,( W ) 2 H, 



i=l 



+ y^TT 3V (H 5 + ^TT 3V H A5 (19) 



where 



= {(#) ll ,(#) 2 ,( W ) ll ,( W ) 2 } 
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With the vacuum alignment in Eq.( 
/ 



it is immediate to derive the mass matrix as follows 
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We note that the Higgs field H 5 (i/45) induces a symmetric (antisymmetric) contribution 

T 2 



to M u . Since the product decompositions (TT)i 1 ~ TiT 2 + T 2 Ti and (TT) 2 ~ 

are both symmetric under the exchange of the two tenplets, the combinations of T 3 T 3 H 45 
and TTif 45 with arbitrary number of flavon fields don't contribute to up quarks masses. 
Therefore the corresponding terms are omitted from the beginning. The mass matrix M u 
is diagonalized by bi-unitary transformation 

V^MuVl = diag(m u ,m c ,m t ) 



The up type quark masses are given by 

2(y C 2^J + ycAV*){y u tiv<t > v x l 'A 2 + y ut 2V v v v> /A 2 ) 2 
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The mixing matrix V£ is 
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We note that there is a mixing of order A 2 between the first and the second family, although 
the (12) and (21) elements of M u vanish at LO. The top quark mass is generated at tree 
level, and the mass hierarchies among the up quarks are reproduced naturally given the 
Vff, and v v of order A 2 A. 



VEVs v x , v v , 



2.3 Down type quarks and charged leptons sector 

The superpotential generating the masses of down quarks and charged lepton is 

9 
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where dots stand for higher dimensional operators. 

o {2) = {x 2 0, x\ <px<f>> vxvi v^V, x0C xvC, </^C, <K 2 } 

0( 4 > = {0 2 X , 0V, 2 C, #X, #C, ^X} (26) 

For the last three terms in Eq. (l25|) . one operator frequently induces several different con- 
tractions, we should take into account all possible independent contractions for each opera- 
tor. Note that the auxiliary symmetry Z 3 x Z 4 imposes different powers of the flavon fields 
for the bottom (tau), strange (muon) and down quark (electron) mass relevant terms. In 
the expansion in powers of 1/A, the bottom (tau) mass is generated at order 1/A, the 
strange (muon) and down quark (electron) masses are generated at order 1/A 2 and 1/A 3 
respectively. As a result, if we only consider the LO operators suppressed by 1/A and 
1/A 2 , the down quark and electron would be massless. 

Recalling the vacuum configuration in Eq. (12|) , we can write down the mass matrix for 
down quarks and charged leptons as follows 

M d =\ 2/ 2 V% 2y d 22 e 2 v IS + y d ' 2 e\ y d 3 e 3 v^ 1 (27) 





yf 2 ^n i 


yf 3 e 3 v^ + 2yf' 3 e 3 v^ 




2?/2V 2w 45 + y^ v z 


y& 3v 5 


2y& 2 m + y!W\ 


yt^\ 


yts ev 5 




yti^5 


-Qy d 2 e 2 v^ + y^e 




-Sy^vis + yi^ 3 ^ 


yt 2 ^ 3 n 


yf 3 e\ - Gyf'ze 3 ^ 


2/23^% 


yl3 £V 5 



M e =\ yf 2 eh^ / -Qy d 22 e 2 v^ + y d 22 e 3 v z y d 2 e% (28) 

where the coefficients yfj = 1,2,3), yf 3 , y 22 and y 31 are linear combinations of the 
leading order coefficients, y 33 coincides with the LO parameter y b up to corrections of 
order e 2 which origins from the operators T 3 FO^ H§. The factor of 3 difference in the 
(13), (22) and (31) elements between M d and M £ is the so-called Georgi-Jarlskog factor [U], 
which is induced by the Higgs H-^. Similar to the up type quarks, the mass matrices M d 
and Mi can be diagonalized by the following transformations 

V d] M d V d = diag(m d , m s , m b ), M e V[ = diag(m e , m„, m T ) (29) 

The mass eigenvalues are given by 

m d ~ \y d xi e^ - y{ 2 y A 2X e\\l \2y d 22 v^) - 2y d 13 y d 22 e 4 v I5 /y d 33 - 4^ 3 y 2 V4r/(^)l 
m s ~ |2^ 2 e 2 W45 + y 2 V%l 

m b - |2/33 £ %I ( 3 °) 
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8 



The diagonalization matrices V£ and V[ are 
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As is shown in Eq. (l30|) and Eq.( l3Tj) . obviously we have 
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(34) 



The well-known bottom-tau unification and the Georgi-Jarlskog relation jH] between the 
down type quark and the charged lepton masses for the second generation are produced in 
the present model. We note that generally the vanishing of the (11) elements of both the 
down quark and charged lepton mass matrices is required, to obtain the Georgi-Jarlskog 
relation for both the first and second generations simultaneously. The mass difference of 
electron and down quark is induced by the Higgs field H-^, acceptable values of the masses 
for electron and down quark can be accomodated due to the Georgi-Jarlskog factor. It 
is well-known that the CKM mixing between the first and the second family is exactly 
described by the Cabibbo angle. In order to satisfy this phenomenological constraint, for 
the parameters yfi an d 2/22 °f order 0(1) we could choose 



J 45 
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Furthermore we can see from Eq. (l32|) and Eq. fl33|) that the mixing angle between the first 
and the second family charged leptons approximately is A c /3. The resulting quark mixing 
matrix is given by 

Vckm = V?V£ (36) 
We can straightforwardly read the CKM matrix elements as follows 
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We note that the CKM elements V us , V c d, V u t and V t d are dominantly determined by the 
mixing in the down type quark sector, V c b and Vu origin from the left handed up quarks 
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mixing. Considering 1*45/% ~ A c , we find that V us and V c d are of order A c , V c b and V ts 
of order A 2 , V u b and V t d are of order A^. The correct pattern of CKM mixing matrix is 
reproduced. 

Due to the non-trivial mixing V[ present in the charged lepton sector, we note that the 
lepton mixing is not the TB mixing, although the light neutrino mass matrix is exactly 
diagonalized by the TB mixing matrix. The lepton mixing matrix (PMNS matrix) is given 
by 

Upmns — V L U v (38) 
Consequently the lepton mixing angles are 
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Taking into account the results for quark mixing shown in Eq. (137j) . we have \V US \ ~ 
1^-— el ~ A c . As a result, the model predicts the deviation of the lepton mixing from 
the TB pattern as follows 
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2.97° 



sin 2 9 



23 



3' 
1 

T 



9 
A| 
36 



A, 



(40) 



The lepton mixing angles are predicted to be in agreement at 3a error with the experi- 
mental data p~H3]. It is remarkable that Eq. poj) belongs to a set of well-known leptonic 
mixing sum rules [42|I43] , and the same results have been obtained in Ref. [28J . 



2.4 High dimensional Weinberg operators 

In the previous section, the neutrinos acquire masses via the see-saw mechanism. It is 
interesting to note that the higher dimensional Weinberg operator could also contribute 
to the neutrino masses directly. In the present model, these effective light neutrino mass 
operators are§ 

w eff = y^ {FFhiXH5H5 + y^ {FF ) 2lp H 5 H 5 + y^(FF) 3lX H, 5 H, 5 + ^(FF) 2( pH 45 H 45 

(41) 

2 Concretely the operator (FF) 3l xH4 5 H 45 denotes yabc{F a ) a {F b )[ ) Xc{H ib )l OL {H 45 )tf , where the Greek 
indices are contracted in the SU(5) space, and the Latin indices are contracted in the ^4 space, the 
coefficient with three S4 indices y a b c is the Clebsch-Gordon coefficient of the S4 group, its value can be 
read directly from the product decomposition rules shown in Appendix A. so that the effective operator is 
invariant under the flavor group S4. The contractions of the remaining operators in Eq.(|41j) can be read 
out similarly. 



10 



Fields 


x° 


T 


6° 

T 


P° 

r 


A 


S 4 


3 2 


2 


3i 


l x 


3i 


z 3 


1 


1 


u 2 


cu 2 




z 4 


1 


1 


-1 


-1 


-1 


U(1)r 


2 


2 


2 


2 


2 



Table 2: Driving fields and their transformation rules under the symmetry group S4 x Z3 x Z4. 



Taking into account the vacuum alignment (x) = v x (l, 1, 1) T and (ip) = 1^(1, 1) T , the 
Weinberg operators in w*$* lead to the following effective light neutrino mass matrix 

/ 2a -a + P -a + p\ 
M e J f = j -a + 2a + /3 -a (42) 
\ -a + p -a 2a + P / 

where 

" = (2^ + 24^ )| 
= (2^1 + 24^^)^ (43) 

The mass matrix M*** is exactly diagonalized by the TB mixing matrix 

U^ B M e J f U TB = diag(mi // , mf, m^) (44) 

As a result, the lepton mixing angles displayed in Eq. (I39[) are not corrected even if the 
Weinberg operators are taken into account. The mass eigenvalues are given by 

eff 



m/ = 3a — P 

h 



m% ff = 2P 



m^f = 3a + P (45) 
Comparing with the see-saw mechanism induced masses in Eq.([n]), we have 



m e ff y2 

(46) 



mi A 2 

It is obvious that the contributions of the Weinberg operator are highly suppressed relative 
to those induced by the see-saw mechanism, so that they can be completely negligible. 



3 Vacuum alignment 

In this section we discuss the minimization of the scalar potential in order to justify the 
vacuum alignment used in the previous section. As usual we introduce a global continuous 
U(1)r symmetry which contain the discrete R— parity as a subgroup. The flavon and 
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GUT Higgs fields are uncharged under U(1)r, the supermultiplets containing the standard 
model matter fields and the adjoint field A carry U(1)r charge +1. Moreover, we include 
additional gauge singlets, the so called driving fields x° ; V 9 ^ 0°; P° an d A with U(1)r 
charge +2. They transform in a non-trivial way under the flavor symmetry S4 x Z3 x Z4, 
as is presented in Tabled Since the driving fields carry +2 unit U(1)r charge, they enter 
linearly into the superpotential. The LO superpotential depending on the driving fields, 
which is invariant under the flavor symmetry, is given by 



«>v = AX O (X^)3 2 +/2X O XC + /3^ O (XX)2 + /4^ O (W) 2 + /5^VC + ^10° 

+<? 2 0°(#)3 1 + gsP^HK + 94P (vv)u + ^iA°(AA) 3l + h 2 A°A£ (47) 

In the SUSY limit, the equations for the minimum of the scalar potential are obtained by 
deriving w v with respect to each component of the driving fields 



8Wy 

M 

dw v 
dw v 

M 

dw v 

dW n , 



fifaiXs - ^2X3) + /2X1C = 
h^PiXx - ^2X2) + /2X3C = 
fi(tpiXa ~ V2X1) + /2X2C = 

fsixl + 2x1x3) + hvl - hfiC = (48) 



This set of equations are satisfied by three types of vacuum alignment 




with the conditions 



v y = T v %-> v <p undetermined (50) 

3/; 



The second is 




(X) = v x I 1 J , (<p) = v v ( _\ ) , (0 = (51) 



with the relations 



2 72/4 + 2/1/2/5 2 f- 2 



2 J2Jt ' "J U«o 2 J * /rr>\ 

v * = — V « = ~W^ (52) 



where is undetermined. The third solution is 




(X) = I I , (<p) = v v ( \ I , (C) = (53) 
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with 

/5 

Vtp = -rV£, undetermined (54) 

J4 

Thus, without assuming any fine-tuning among the parameters = 1 — 5), the VEVs 
v x and v 9 are expected to be of the same order of magnitude for the three cases, 

v x ~ v v (55) 

Only the first alignment can produce the results in the previous sections, we need of some 
soft masses in order to discriminate it as the lowest minimum of the scalar potential, 
since the values of the scalar potential for the three solutions are exactly the same in the 
SUSY limit. It is well-known that the soft mass usually is of order TeV, consequently 
the difference of the scalar potential for different vacuum solutions is marginal comparing 
with the flavon VEVs. As has been shown in the previous section, at LO the S4 flavor 
symmetry is broken by the VEV of \ and ip in the neutrino sector. The flavon fields 0, 
r), A and £ are involved in generating the quark and charged lepton masses, their vacuum 
configurations are determined by 



dw v 



<90° 
dw v 

M 

dw 



201 (01 - 0203) + fl'2('7l02 + ??203) = 
2#l(02 - 010 3 ) + #2(771^1 + ??202) = 
2#l(03 - 0102) + #2(?7l03 + ??20l) = 

- g 3 (0? + 20 2 3 ) + 2g mm = (56) 



dp 

dw. 



- 2h l (A 2 1 - A 2 A 3 ) + h 2 A^ = 
2h 1 (A 2 2 -A 1 A 3 ) + h 2 A 3 ^ = 



2^ 1 (A^-A 1 A 2 ) + /i 2 A 2 e = (57) 



The equations in Eq. (l56|) lead to two different vacuum configurations, the first is 



(0) = v, , (rj) = U ) (5 




with 



The second is 



Vcf, = —-^—v v , v v undetermined (59) 



2g 
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with = and v v undetermined. Obviously and v v are expected to be of the 

same order for both solutions. As before, we select the first vacuum configuration. The 
equations in Eq. fl57|) admit three un-equivalent solutions, the first is 



(A) = , (0 = v e (61) 




with t>A = ~y^ v £ an d undetermined. In a similar way, the VEVs v& and should be 
of the same order. The second solution is 




(0 = (62) 
with v& undetermined. The third vacuum configuration is 

(A) = ( ) , (0 = v< (63) 




and remaining undetermined. The first vacuum solution is chosen in the present model. 
Acting on the vacuum configurations of Eq.( l4*9l) . Eq.( l58l) and Eq.( l6TT) with the elements 
of the flavor symmetry group S4, we can generate other minima of the scalar potential. 
These new minima are physically equivalent to the original set, and they all lead to the 
same physics, i.e. fermion masses and flavor mixings. As a result, without loss of gener- 
ality we can analyze the model by choosing exactly the vacuum in Eqs.( )49ll58l[6TT) as local 
minimum, and the different scenarios are related by field redefinitions. Since no superpo- 
tential couplings of positive mass dimension are involved in the flavon superpotential, the 
trivial solution with all flavon VEVs vanishing can not be excluded. However, by taking 
into account the contribution of the soft mass terms, we can discriminate the configuration 
in Eqs. ()49ll58|l6"Tj) as the lowest minimum of the scalar potential. Regarding the size of 
the flavon VEVs, we have the relations v x ~ v v , ~ v v and i>a ~ v^, as have been 
demonstrated above 0. Furthermore, The magnitudes of the flavon VEVs are determined 
by the patterns of fermion mass hierarchy and mixing. From the LO predictions presented 
in the previous sections, we find that in order to reproduce the correct patterns of fermion 
masses and flavor mixings, a common order of magnitude for the VEVs scaled by the 
cutoff A is expected 

Vx V ± V _H Vi ,2 / fi4 N 

A~A~A~A~A~A C { } 



We consider the soft terms involving \ ancl which is generally written as m£|x| +m^\ip\ +m^x + 
m^ip 2 . By choosing rr? x , m^, and rh 2 v < 0, the vacuum shown in Eqs. (|49l58l61l) are more stable than 
the vanishing VEVs configuration. 

In the absence of specific dynamical tricks (see Ref. |Sj for a model in which such a trick is imple- 
mented), the uncorrelated VEVs naturally have the same order of magnitude, this is consistent with 
the the constraints from the measured mass hierarchies and flavor mixing. Moreover, we note that the 
correlation of scales of more flavon VEVs can be achieved by adding further driving fields, whereas this 
procedure would introduce more fields and free parameters. 
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Moreover, we will show in the following that the successful LO predictions are not destroyed 
by the subleading corrections when the vacuum alignment is chosen as has been stated 
above. Similar to the fact that the F-terms of the driving fields are the origin of the 
alignment of the flavon VEVs, we can derive the vacuum structure of the driving fields 
from the F-terms of the flavon fields. As all terms in the flavon superpotential are linear 
in the driving fields, the configuration in which all these fields have vanishing VEVs is in 
any case a solution. Moreover, we have checked that this is unique vacuum configuration 
of driving fields in our model. 



4 Subleading corrections 

At the next level of expansion in 1/A, the superpotentials w u , w u , Wd and w v are 
corrected by higher dimensional operators whose contributions are suppressed by at least 
one power of 1/A. The corrections to w v result in small deviations from the LO vacuum 
alignment thus affect the results for fermion masses and mixings. In addition the fermion 
mass and mixing matrices are corrected by the subleading operators in w u , w u and Wd- 
In the following we shall first present the analysis for the subleading corrections to the 
vacuum alignment, then move to the corrections to fermion mass matrices. 



4.1 Corrections to the vacuum alignment 

We detail the discussion of this issue in Appendix B, here we only present the results. 
The vacuum configuration is shifted into 

(X) = I v x + Sv x2 | , <*>=( .Xj, <0 = *»c 




(0) = v 4 , + Svja , {77) 



8v vl 



(A) = 5v A2 , (0= vs (65) 



where v v , v v and remain undetermined. Since we typically have VEV/A ~ \ 2 C at LO, 
we expect the order of magnitude of the shifts as follows 

&W u X~ A c> Sv x3/ V X~ X li HsK~ A c> Sv </ v X~ X l 

Svai/va ~ A^, 5va2/va ~ A c , 5v^/va ~ A^ (66) 

From Appendix B, we can see that the subleading operators linear in x° an d y^ are 
suppressed by A, while the subleading operators linear in 0°, p° and A are suppressed by 
A 2 , due to the constraint of the flavor symmetry S4 x Z 3 x Z4. Consequently the shifts 
are of order \ 2 C or A^ with respect to the LO VEVs, as is displayed in Eq.( 
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4.2 Corrections to fermion masses and mixings 

The fermion masses and mixings receive corrections from both the shifted VEVs and 
the higher order terms in the superpotentials w u , w u and Wd- We can obtain the mass 
matrix by inserting the modified VEV into the LO operators plus the contribution of the 
higher dimensional operators evaluated with the LO VEVs. For convenience, we denote 
X, v 9 an d C with <& u , A and £ with $fi, <p and i] with $y 2 in the following. 

4.2.1 Corrections to up quark sector 

First we discuss the corrections to the up type quark mass matrix coming from the 
modified vacuum alignment. Plugging the shifted vacuum configuration shown in Eq. (l65j) 
into the LO superpotential in Eq. (]19|) . we find that the (12) element receives a correc- 
tion of order e A from the operators (TT)i 1 ((j)(())i 1 H 5 and (TT) ll (r)T])i 1 H 5 . The operator 
(TT) 2 (00)2-f^5 induces corrections of order e A to both (11) and (22) elements. The cor- 
rections to the (13) element due to modified VEVs are of order e 3 , they arise from the 
contractions TT 3 (0x) 2^/5, TT 3 (rjip) 2 H 5 and TT 3 r]H^. Meanwhile, TT 3 (0x) 2^/5 introduces 
correction of order e 3 to the (23) element. 

Then we come to discuss the corrections caused by the higher dimensional operators 
in the matter superpotential w u . We note that all corrections to the (33) element can be 
absorbed into the coupling y t of the LO operator T 3 T 3 H 5 , and any operator comprising 
the superfields T 3 T 3 H^ (TTif 45 ) and arbitrary number of flavon fields gives a vanishing 
contribution to the up quark mass matrix. Due to the auxiliary symmetry Z 3 x Z 4 , the 
subleading corrections to (11), (12), (21) and (22) elements arise at order 1/A 4 , they come 
from the following contraction 

^IT$A$/2V6 ( 67 ) 

Consequently the corrections to the = 1,2) element from the high dimensional 

operators are of the same order as those from the shifted vacuum. In the same way, we 
find the (13) and (23) elements are corrected by the following contractions 

^TT 3 $ fl ® f2 $ f2 H 5 , ^TT 3 <$> v <$> v $ f2 H i5 (68) 

Substituting the LO VEVs into the above contractions, we notice that the corrections to 
the (23) and (32) elements originate from the latter contraction, whereas both operators 
contribute to (13) and (31) elements. In short summary, the up type quark mass matrix are 
corrected by both the deviations from the LO VEV alignment and the higher dimensional 
operators allowed by the flavor symmetry. We can parameterize the up quark mass matrix 
as 

/ Sy^e 4 v 5 8y& 4 v 5 ±y& 2 v 5 + 8y£e 3 v 45 \ 

M u =\ 8y? 2 e 4 v 5/ Sy^e 2 v 5 ^ 8^45 + 4<y&e 2 v 5 (69) 

where y^ (i, j = 1, 2, 3) are complex numbers with absolute value of order one, and they are 
linear combinations of the leading and subleading coefficients. We note that y 33 coincides 
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with the LO parameter y t up to higher order corrections which are due to two flavons and 
three flavons insertions in the operator T 3 T 3 H 5 . Similarly the parameters ?/" 3 , y 22 , 2/23 anci 
2/23 are determined by the LO couplings in Eq. (1191) up to small corrections of relative order 
e or e 2 . The mass matrix M u in Eq.( |69l) leads to the up type quark masses 



m,, 



V22 V 5 



(2/23, 



]y ' lV5 2 vh V5 + y^(y3s) 2 vi + yUyh) 2 v 2 



45 

2 „,2 



?/33 V 5 

m t ~ |8t/ 3 > 5 | (70) 
The correct hierarchies among the up quark masses are obtained. 

4.2.2 Corrections to down quark and charged lepton sector 

Plugging the shifted vacuum of <fi into the LO operator T 3 F<pH-^ leads to corrections 
to the (13), (23) and (33) elements of M d of order e 3 %, this amounts to a rescaling of 
the parameters yf 3 , y 23 and y 33 . If the non-zero shifts Svai,2,3 are taken into account, the 
LO operators TFAAHj^ and TFA^Hj^ introduces corrections of order e a vj§ to the (22), 
(31), (32), (11), (12) and (21) elements which also receive corrections of order e 4 v^ from 
the operators TFO^'H-^. At LO the superpotential Wd is expanded to 1/A 3 , it is corrected 
by the following subleading operators 

1 1 1 

—T 3 F^ fl ^ f2 ^ f2 H^, -^TFQf&ftQfzQfiHjr, —TF$ fl $ fl $ fl $ fl Hjr, 
^TF^^ fl ^ fl H^ (71) 

With the LO VEVs, the above high dimensional operators lead to corrections of order 
e 4 v^ or £ 4 t>45 in each entry of Md- Note that the subleading terms involving the fields 
combination T 3 FH^ arise at order 1/A 5 with the insertion of five flavon fields. Therefore 
we conclude that the NLO corrections to the down quark and charged lepton mass matrices 
can be reabsorbed into a redefinition of the LO parameters whose order of magnitudes are 
not changed. As a result, we can parameterize the down quark and charged lepton mass 
matrices in the same way as Eq.f l27|) and Eq. fTSHl) . and the parameters yfAi,j = 1,2,3) 
are still used to avoid introducing extra unnecessary parameters. However, we should 
keep in mind that the value of yfj is different from the corresponding LO one due to the 
subleading corrections of relative order e or e 2 . Since the down quark and charged lepton 
mass matrices remain the same form, the down quark and charged lepton masses are still 
given by Eq. (!30l and Eq. (l3T|) respectively. Starting from the quark mass matrices M d 
in Eq.f l27p and M u in Eq. f !69j) . we can straightforwardly find the CKM mixing matrix as 
follows 



4d - Vcs ~ Vtb 



Vil V 5 , 1 J/l 3 J/23 u 5V45 



Vus ~ 2<4 v^" ' 2 y^ 3 v 2 + (y u 23 ) 2 v 2 , 
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V ub — Z d £ 



vii 2/13 , i y u M) 2 vl 5 



yis %33 2 yUy^l + vh(vh) 2 < 



5 2 



2/33 ^5 



,<2 



4, ~ -2^£ 
2/33 % 



2/212/23 ^45 n 2/fi , 2/is n 



22/222/33 ^5 ^45 2/33 ^33 



e 2 (72) 



We see that the CKM parameters are determined by the LO results up to small corrections, 
which is absorbed in the redefinition of the parameters. The successful LO predictions for 
the order of magnitudes of the CKM matrix elements are not spoiled by the subleading 
corrections. 

4.2.3 Corrections to neutrino sector 

The superpotential w u in Eq. (j3J) is corrected by terms with more flavons insertion. The 
VEV shifts in the LO operator do not affect the Dirac mass, the NLO corrections to the 
neutrino Dirac couplings are 

^l(FA) 3l xH i5 + y^(FA) 2( pH 45 (73) 

Using the alignment of x an d <p given in Eq.( l49j) . the NLO corrections to the Dirac mass 
matrices read 

^Vv\± ~Vv\-K + 2/^x -2/^ix + 2/^x 



5M° = ~ V45 1 -y vl % + yv2 % 2y„i v it + y»i\ -y*i 

-Vvi x + 2/^2 x -2/^1 x 2 Vvif + y»2 



\/1 5 

6M° = ~SM^ (74) 

We note that the subleading corrections 5M^ 3 and 5M® are still compatible with the TB 
mixing. The Majorana mass matrices of p 3 and po are modified by the terms 

X l A 2 5 X + A 2 A 2 ^ + X -j-A\H 2A + ^A 2 <pH M + ^(AA^A^ + ^f(AA) 2 (<pA) 2 
+^(A4) 3l (0A) 3l + ^-(AA) 3l ^+^(AA) 3l ( V A) 3l + ^(AA) 2 t£ (75) 

where 5\ and Sip denote the shifted vacua of the flavons x an d V 9 - The operators A 2 xH 2 ^ 
and A 2 tpH 24: lead to mass splitting between the triplet p 3 and the singlet po, and the 
resulting contributions to the Majorana mass matrices have the same structure as the LO 
predictions. Taking into account the possibility of absorbing the corrections partly into the 
LO results, the remaining operators give rise to three independent additional contributions 
to the Majorana mass matrices. As a result, the next to leading order corrections to 
and can be parameterized as 



—6xa€V x 3xa£V x — 3xB^V tp 3xacv x — 3xB^v tp 

SM^ = | 3xacv x — 3xbcv v —Gxa£v x — 3xbcv v 3xa^v x 

3xacv x — 3xb€V v 3xacv x —6xacv x — 3xb£v v 
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SM 



M 
PO 



(0 -x c x D - x E 

-x c x D + 2x E 
x D - x E 2x c 

—2xa£v x XA^v x —x B ev ip 
xa£v x — x B ev^ —2xa£v x — x B ev v 
xa£v x — x B ev^ xa£v x 

-x c x D - x E 

—x c x D + 2x E 
xd — %e 2xc 



x 
A 



(76) 



x A ev x - xsevp 

xa^v x 
-2xa£v x — x B ev v 



X 

A 



(77) 



where e = ^t=^IE The first term in both Eg. (176 P and Eq. fl77|) represents the contribution 
of ^A 2 xH24 and ^-A 2 ipH2A, the second term denotes the effects of the modified vacuum 



configuration and the sub leading terms with the form AA& /i$/ 2 , which leads to deviations 
from the TB mixing pattern in the neutrino sector. With the LO contributions shown in 
Eqs. fHfo'j) and the NLO corrections in Eqs. (l74|76|r77|l . the light neutrino mass matrix can 
be obtained straightforwardly via the see-saw formula. To first order in small parameters 
v x /A, v^/A and e, we find that the light neutrino masses are given by 



mi 



m 2 



m 3 



9(-3x A v x + x B v v )ylv 2 ^ (2x c - x D + 2x E )v 2 x ylvl 



5(-3Aiw x + A 2 f ¥ ,) W(-3Xiv x + \ 2 v^y 



-e + 



h(-3\ lVx + A 2 ^) 2 A 



vl4 



9x B y 2 u vf 



20\ 2 2 v v 



< 5 e x D v\ylv\ 



10A|^A 



9{3x A v x + x B v v )ylvl (2x c + x D 



2x E )v x ylvl 



5(3A lUx + A 2 ^) 2 A 



5(3Ai^ x + A 2 i^) 10(3Aif x + A 2 ^) 2 
The lepton mixing angles are modified as 

x D v 2 x (3\\v* x + A*^) + (3A^ X - \ 2 v v )S* D (v 



sin 9 



13 



f 3 



-e - 



*\2 
X> 



22 U 45 



(x c - x E )v 2 (3Xlv* + \*v*) + 2\ 2 v v (x* c - x* E ){v*) 



y/2 [| - 3Xn> x + A 2 tvl 2 - \3Xiv x + A 2 ^| 2 ]A 

*\2 



\/2[4|A 2 ^| 2 - |3A^ X + A 2 ^| 2 ]A 



sin 2 9 12 



1 

~ 3 



i8 [ y 2 W 

.2, 



2/22 V 45 



3[| -3A^ X + A 2 ^| 2 -4|A 2 ^| 2 ]A 



sin 2 9 



23 



x[(x c + x E )v x (3\lv x 



1 
2 



(^ + ^)K) 2 (3A lWx + A 2 t V )] 



3x D v x \l 


v x 


2 + 3#bv* Ai 




2 


2[|3A^ X + A 2 ^| 2 




- 3Aiv x - 


1- A 2 v v ,| 2 ]A 



2[4|A 2 ^| 2 - 1 3 



x [(x c - x s K(A>* + A;<) + (x c - ^)K) 2 (A^ X + X 2 v v )] 



A 2 ^| 2 ]A 
(79) 



We note that the subleading terms in the neutrino sector induce corrections of order e to all 
three mixing angles, and the deviation of 9 2 % from its TB value is mainly determined by the 
NLO contributions in the neutrino sector. The lepton mixing angles are still compatible 
with the current experimental data. In particular, the reactor angle #13 is within the reach 
of next generation neutrino oscillation experiments. 



5 The SU(5) GUT symmetry is broken down to the standard model symmetry by the VEV of f/24, i.e. 
(tf 24 ) - v-u/VSO diag(2, 2, 2, -3, -3). 
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Figure 2: The allowed region of sin 2 6\2 and sin 2 023 against sin 2 #13 for both normal hierarchy and inverted 
hierarchy neutrino mass spectrum 



5 Phenomenological consequences 

In this section, we shall present the predictions for some phenomenologically interesting 
observables in our model. Here we are particularly interested in the neutrino sector, both 
the LO and the NLO contributions are taken into account in the following. 

5.1 Lepton mixing angles 

The analytic expressions for the mixing angles are shown in Eq. (J79l) , which allows us to 
estimate the deviations from the TB mixing pattern qualitatively. In order to see clearly 
the allowed ranges of the mixing parameters, as well as for cross-checking the reliability 
of the analytical estimates, we shall perform a numerical analysis. All the involved LO 
and NLO coefficients are taken to random complex number with absolute value in the 
interval [1/3, 3], the expansion parameters e and e have been fixed at the representative 
value of 0.04 and 0.001 respectively §, and the VEV ratio ^45/% is set to the indicative 
value 0.22. Furthermore, we require the oscillation parameters Am^, Am^, sin 2 # 12 , 
sin 2 #23 and sin 2 #13 to lie in their 3cr interval. The allowed regions of sin 2 612 — sin 2 #13 
and sin 2 #23 — sin 2 #13 for both normal hierarchy and inverted hierarchy are showed in Fig. 
[21 It is obvious that rather small #13 is favored for both NH and IH spectrum, which is 

6 For other indicative values of the small parameters e and e, the resulting numerical results don't 
change qualitatively. 
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consistent with our theoretical analysis. 

5.2 Neutrinoless double beta decay 

The discovery of neutrinoless double beta decay 0^2/3 is very important because it 
could directly establish lepton number violation and the Majorana nature of neutrino. The 
0z/2/3 decay amplitude is proportional to the effective mass \m ee \, which is (11) element of 
the neutrino mass matrix in the basis where the charged lepton mass matrix is real and 
diagonal. The allowed region for m ee is displayed in Fig. El The horizontal lines denote 
the future sensitivity of some Oz/2/3 decay experiments, which are 15 meV and 20 meV 
respectively of CUORE [15], Majorana [I6]/GERDA III [U] experiments. 
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Figure 3: Scatter plot of the effective mass \m ee \ with respect to the lightest neutrino mass. The blue 
corresponds to the normal hierarchy neutrino spectrum and the red to the inverted hierarchy spectrum. 

We see that \m ee \ is predicted be below the present bound from the Heidelberg- Moscow 
experiment [4"4"] . For the NH spectrum, \m ee \ can be so small to be close to zero, while 
the scatter plot indicates a lower bound for \m ee \ of about 14 meV in the case of IH 
spectrum. It is quite close to the future experimental sensitivity so that 0z/2/3 decay should 
be observed by future experiments for IH spectrum. We note that a partial cancellation 
between the LO and NLO contributions takes place, so that the lightest neutrino mass can 
be very small of order 10 -4 eV in contrast with the LO constraints shown in Eq. fll7|) . This 
however requires an additional fine tuning of the parameters which has been reproduced 
in our numerical analysis only partially and by very few points. 

In FigHl we plot the effective mass mp = [^2 k \(UpMNs)ek\ 2, m'l] 1 ^ 2 in /3 decay experi- 
ments, which could measure the non-zero neutrino masses. We clearly see that the effective 
mass 771/3 is predicted to be below the prospective sensitivity 0.2 eV of the KATRIN ex- 
periment for both NH and IH neutrino mass spectrum [4~8] . 
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Figure 4: mp as a function of the lightest neutrino mass. The blue corresponds to the normal hierarchy 
neutrino spectrum and the red to the inverted hierarchy spectrum. 

5.3 Sum of neutrino masses 

The prediction for the sum of neutrino mass is presented in Fig. [51 The horizontal line 
is the cosmological bound at 0.19 eV, which is obtained by combining the data from the 
Cosmic Microwave Background (CMB) anisotropy (from WMAP 5y |49j . Arcminute Cos- 
mology Bolometer Array Receiver (ACBAR) [50], Very Small Array (VSA) [51J, Cosmic 
Background Imager (CBI) [52] and BOOMERANG [53] experiments) plus the large-scale 
structure (LSS) information on galaxy clustering (from the Luminous Red Galaxies Sloan 
Digital Sky Survey (SDSS) [51]) plus the Hubble Space Telescope (HST) plus the lumi- 
nosity distance SN-Ia data of [55] plus the BAO data from [56] and finally plus the small 
scale primordial spectrum from Lyman-alpha (Lya) forest clouds [57]. We see that our 
model predicts J2k m k too similar for both hierarchies to be distinguished using the current 
cosmological information on the sum of the neutrino masses. 

The present model has rich phenomenological implications, we only study few interest- 
ing observables here. Especially the predictions for lepton flavor violation branching ratios 
and leptogenesis deserve to be studied carefully, which are important to test the model 
and distinguish this model from other discrete flavor symmetry models. These topics will 
be discussed in future work [58] . 

6 Conclusion 

In this work, we have built a SUSY SU (5) model based on the flavor symmetry S4 x 
Z 3 x Z4. Three generations of adjoint matter superfields are introduced, and they are 
assigned to transform as 3i of S4 . The neutrino masses are generated via the combination 
of type I and type III see-saw mechanism in the model. To describe quarks, we make use 
of the doublet representation of 5*4, we accommodate the first two generations of tenplets 
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Figure 5: The sum of neutrino masses J2k mk as a function of the lightest neutrino mass. The blue 
corresponds to the normal hierarchy neutrino spectrum and the red to the inverted hierarchy spectrum. 



10 in doublet under S±. In particular, the first generation 10 is assigned to the second 
component of the doublet, and the second generation as the first component, whereas 
the third generation of 10 is kept invariant. The observed mass hierarchies of quarks 
are reproduced naturally via the spontaneous breaking of the flavor symmetry without 
invoking the Froggatt-Nielsen mechanism. In order to generate the CKM mixing between 
the first two generations, we require a moderate fine tuning v-^/v-^ ~ A c , then the model 
generates the observed pattern of the CKM mixing matrix. 

In the neutrino sector, the flavor symmetry S4 is broken down to the Klein four sub- 
group by the VEV of the flavon fields x an d <p at LO. The resulting light neutrino mass 
matrix is exactly diagonalized by the TB mixing matrix, and the neutrino mass spec- 
trum can be normal hierarchy or inverted hierarchy. There are only three independent 
parameters in the neutrino sector at LO, the model is rather predictive so that the lightest 
neutrino mass is constrained to be larger than 0.011 eV and 0.028 eV for normal hierarchy 
and inverted hierarchy respectively. The mixing of the left-handed charged leptons results 
in corrections to the TB mixing pattern, which is described in terms of a well-known lepton 
mixing sum rule, and the reactor mixing angle is predicted to be close to three degrees. 

The subleading corrections to the flavon alignment and the fermion mass matrices have 
been analyzed carefully. We show that the successful LO predictions for the pattern of 
quark masses and CKM mixing angles are not spoiled by the subleading contributions, and 
all the three lepton mixing angles receive corrections of order A^. The phenomenological 
implications of the model are investigated in details, we find that future neutrinoless double 
beta decay experiment with high precision is an important test of the model, it allow us 
to distinguish the NH spectrum from the IH one. Finally we note that the predictions 
presented in the work are valid just below the GUT scale. To determine the fermion masses 
and mixings and phenomenologically interesting observables at the electroweak scale, we 
should study the renormalization group running carefully. These issues will be studied 
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elsewhere in future [58J. 
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Appendix A: The group S4 and its representation 

The group £4 is the group of the permutations of four objects, it has 4! = 24 elements. 
Let a generic permutation be denoted by (1, 2, 3, 4) — > (m, ri2, n 3 , 77.4) = {ni^nsn^) . £4 can 
be generated by the two basic permutations £ and T where S = (2341) and T = (2314). 
We can check that £ 4 = T 3 = 1, ST 2 S = T. £ 4 has five conjugate classes as follows 

d : 1 

C 2 : STS 2 = (2134), TSTS 2 = (3214), ST 2 = (4231), S 2 TS = (1324), TST = (1432), 

T 2 S = (1243) 
C 3 : TS 2 T 2 = (2143), £ 2 = (3412), T 2 S 2 T = (4321) 

C 4 : T = (2314), T 2 = (3124), T 2 £ 2 = (2431), S 2 T = (4132), S 2 TS 2 = (3241), £T£ = (4213), 

S 2 T 2 = (1342), T5 2 = (1423) 
C 5 : S = (2341), T 2 5T = (2413), ST = (3421), TS = (3142), TST 2 = (4312), S 3 = (4123) 



The structure of the group S4 is rather rich, it has thirty proper subgroups of orders 1, 
2, 3, 4, 6, 8, 12 or 24. Concretely, the details about the subgroups of S 4 can be found in 
Ref . [23] . For a finite group the number of irreducible representation is equal to the number 
of conjugate class. Consequently the S4 group have five irreducible representation: li, 1 2 , 
2, 3i and 32, which are all real. Concretely the representation matrix can be chosen as 




T 




for 


li 


for 


1 2 


'■) 


for 2 


N 







| for 


OJ J 




N 







for 


w J 





S=±\ 2u 2u 2 -1 , T 



5 = ^ 2w 2w 2 -1 , T 



The characters, i.e. the trace of the representation matrix, are given in the character table 
(please see Table [3]). From the character table of the £4 group, we can straightforwardly 
obtain the multiplication rules between the various representations 

li (g) lj = l((i+j) mod 2)+l, lj ® 2 = 2, 1* <S> 3j = 3((i +J -) mod 2)+l 

2 (8) 2 = li © 1 2 © 2, 2 © 3^ = 3i © 3 2 , 3i © 3i = li © 2 © 3i © 3 2 , 
3 X ©3 2 = 1 2 ©2©3i©3 2 , with ij = 1,2 (80) 

Starting from the explicit matrix representation, we get the product decomposition rules 
of the £4 group. In the following we use aj to denote the elements of the first representation 
of the product and to indicate those of the second representation. 
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Classes 




^2 


Co 


Ca 






1 

± 


1 

± 


1 


1 

± 


1 

± 


12 


1 


-1 


1 


1 


-1 


2 


2 





2 


-1 





3i 


3 


1 


-1 





-1 


3 2 


3 


-1 


-1 





1 



Tabic 3: Character table of the £4 group. 



1 2 © 1 2 = li 



1 2 ® 2 = 2 



1 2 (8) 3i = 3 2 



a/3 (81) 



a Pi 
-a/3 2 



1 2 (8) 3 2 = 3i 



a/3; 



3 



2©2 = l 1 ffil 2 ©2 



«2/32 

ai/3i 

3i <g> 3i = 3 2 ® 3 2 = li © 2 © 3i © 3 2 

li ~ ai/3i + a 2 /3 3 + a 3 /3 2 

a 2 /3 2 + ai/3 3 + a 3 /3i 
«3^3 + «i^ 2 + a 2 ^i 



52) 



a/3i 
a/3 2 
a/3 3 

a/3! 

3i ~ I a/3 2 I (84) 



11 ~ ai/3 2 + a 2 /3x (85) 

1 2 ~ ai/3 2 - a 2 /3i (86) 



57) 



59) 



2ai/3i - a 2 /3 3 - a 3 /3 2 
2a 3 /3 3 - ai/3 2 - a 2 Pi 
2a 2 /3 2 - ai/3 3 - a 3 p 1 

«2/3 3 - a 3 /3 2 
aiP 2 - a 2 Pi 
a 3 pi - ai/3 3 
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• 3i <g> 3 2 = 1 2 © 2 © 3i © 3 2 



2 © 3i = 3i © 3 2 



2 © 3 2 = 3i © 3 2 




- aift + a 2 /3 3 + a 3 /3 2 (92) 

a 2 /3 2 + ai/?3 + a 3 £i \ , 
-a 3 /3 3 - ai/3 2 - a 2 /?i / 

«i/3 2 - a 2 A 
a 3 /3i - «i/3 3 

2«i/3i - a 2 f3 3 - 

2a 3 /3 3 - ax/32 - a 2 /3i | (95) 
2a 2 /3 2 - ai/3 3 - 



ai/3 2 + a 2 /3 3 
ai/3 3 + a 2 /?i 
+ a 2 /3 2 

ai/3 2 - a 2 /3 3 
«iA - a 2 /3 2 



ai/3 2 - a 2 /3 3 

3i ~ I ax/33 - a 2 A I (98) 



2 

«i/3 2 + a 2 /3 3 
«i/3 3 + «A 
cci/?! + a 2 /3 2 

We note that the multiplication rules presented above are in accordance with the results 
in Ref. [231. 



Appendix B: Vacuum alignment beyond the leading 
order 

In this appendix we discuss the subleading terms of the driving superpotential w v and 
the resulting corrections to the LO vacuum alignment. At the next level of approxima- 
tion the LO driving superpotential is corrected by operators of higher dimension whose 
contributions are suppressed by at least one power of 1/A. As a result, the superpotential 
depending on the driving fields becomes 

w v =w Q v + Sw v (100) 
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h <</>>=( w * J. (0 = o, with< = -^-< 




The leading order term wJJ reads 

< = /iX°(X^3 2 + f2X°X( + /V(xx) 2 + /VW) 2 + /s^VC + <7i0°( 

+</ 2 0°(#) 3l + <feP°(#)ii + + ^iA°(AA) 3l + /i 2 A°A£ (101) 

In the SUSY limit, u>° leads to the following vacuum configuration 

(0) = | ^ | , fa) = ( J J ' With u * = 

(A) = | | , (0= Uf, with WA = - (102) 

The correction term consists of the most general subleading operators linear in the 
driving fields, and they should be invariant under the flavor symmetry 5*4 x Z 3 x Z 4 . 

1 15 1 10 1 39 i 18 1 28 

^. = xE *<°? s + x E »<°r° + ^ £ «°f + ^ £ + ^ E ( <°. A " ( 103 > 

i=l i=l i=l i=l i=l 

where ki, Wi, Sj, rj and tj are order one coefficients, their specific values are not determined 
by the flavor symmetry {Of , Of , Of , Of , O t A °} denote the complete set of subleading 
contractions invariant under £4 x Z 3 x Z4. 

Of = (X°X) 2 (0A) 2 , of = (x°x) 3l (0A) 3l , of = (X°X)3 2 (0A) 32 , 
of = (xVk(0A) 3l , of = (xV) 32 (0A) 32 , of = (x°xk0£, 
of = of = (x°x) 3l (^A) 3l , < = ( X °x) 32 (r/A) 32 , 

< = (xV) 3l faA) 3l , ot = (xV) 32 faA) 32 , of = (x°x) 2 ^ 

< = X°C(0A) 3l , oS = X°C0£, < = X °C(^A) 3l (104) 



of = (<Ak(0A) 3l , Of = (^°x) 32 (0A) 32 , Of = ( ¥> V)i 1 (^A) ll , 

Of = (<A) 2 (0A) 2 , Of = (<AM£, Of = (^ x) 3l (r/A) 3l , 

of = (A)s a (i?A) aa , Of = (ifPiphtf, Of = ^°C(0A) 2 

Of = ^ (105) 



of = ^(xxkO^k, Of = 0°((xx) 2 (#) 3l ) 3l , Of = 0°(xxk(#)ii, 

of = 0°((xx) 3l (00) 2 ) 3l , of = 0°((xxk(00) 3l ) 3l , of = ^(xxkfa^k, 
of = <A(xx) 2 (#kk, of = 0°((XX)2(#) 32 ) 31 , of = 0°((xx) 3l (#) 3l ) 3 
Of = 0°((xxk(#kk, Of = ^(xxMwk, Of = 0°((XX)3 1 (W) 2 )3 1 , 
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of = ^W)ix(Msx, OH = o ((w) 2 (00) 3l ) 3l , Of, = ^KM-k, 

Of 6 ° = 0°(( W ) 2 (#) 3l ) 3l , Of 7 ° = 0°((W)2(#)3 2 ) 31 , < = 0°^X) 3l (#)!,, 

< = o ((^x) 3l (#) 2 ) 3l , o£ = 0°(^x) 3l (#) 3l ) 3l , < = 0°((^) 32 (<W>) 2 ) 3l , 

< = 0°((^X) 32 (00) 3l ) 3l , 0£ = 0°((^X) 3l (#) 3l ) 3l , OH = 0°((^) 3l (#) 32 ) 3l , 

o& = 0°((^x) 32 (#) 3l ) 3l , < = 0°((^X) 32 (#) 32 ) 31 , Og = 0°(^x) 3l (w)i 1 , 

< = 0°((^X) 31 (W) 2 ) 31 , o£ = 0°((^X) 32 (W) 2 ) 31 , Of = 0°(x(#) 2 ) 32 C 

< = 0°(x(00) 3l ) 32 C, pg = 0°(x(#) 3l ) 32 C, oil = 0°(x(#) 32 ) 32 C, 

ci° = o (x(w) 2 ) 32 c, = 0°(^(#) 3l ) 32 c, o£ = 0°(^(#) 3l ) 32 c 

C# = 0°^(#)3 2 ) 32 C, OH = 0°(#) 3l C 2 , eg = 0°(#) 3l C 2 (106) 

< = AxxM^k, c 2 p ° = P°(XX) 2 (00) 2 , of = p°(xx) 3i (#) 3i , 
of = p°(xx) 3i (#) 3i , eg" = AxxkCwk, 6 p ° = p°(xx) 2 (w) 2 , 

O? = pOWuMu, Of = p°( W ) 2 (#) 2 , < = Awlu^u, 

< = p°(w) 2 (w) 2 , og = pVxk(#) 3l , ^ = p°(<px) 3i (#) 3i , 
p 3 = pVx) 32 (^) 32 , < = P°xC(#) 32 , o£ = pVC(#) 2 , 

p 6 = pVC(w) 2 , C^ = AW)u, o& = (107) 

Of° = AW^AAk, Of = A°(( xx ) 2 (AA) 3l ) 3l , A ° = A°(xx)3 1 (AA) ll , 

Of = A (Oa) 3l (AA) 2 ) 3l , of = A°(( xx ) 3l (AA) 3l ) 3l , of = A^xxkAC, 

O 7 A °=A (( xx ) 2 A) 3i e, Of = A°((xx) 3l A) 3i e, Of = A°(xxk^, 

A ° = A°((^x)3i(AA)i 1 , 0^ = A ((^ x ) 3l (AA) 2 ) 3l , A ° = A ((^ x ) 3l (AA) 3l ) 3l , 

Of = A°((^ x ) 32 (AA) 2 ) 3l , A ° = A°((^x) 32 (AA) 3l ) 3l , Of = A°((^) 3l A) 3i e, 

Cf 6 ° = A°((^x) 32 A) 3i e, Of = A°^xhM, Of = A°((f(p) ll (AA) 3l , 

Of = A°(( w ) 2 (AA) 3l ) 3l! Of = A'i^AC, Of = A°(( w ) 2 A) 3i e, 

Of = A°( X (AA) 2 ) 32 C, 2 A 3 ° = A°( x (AA) 3l ) 32 C, Of = A°( X A) 32 £C, 

O 2 A 5 =A°(^(AA) 3l ) 32 C, Of = A°(^A) 32 eC, Of = A°(AA) 3l C 2 , 

Of = A°A£C 2 (108) 



The subleading contribution 5w v induces shifts in the LO VEVs shown above, then the 
new vacuum configuration can be parameterized as 




v x + 6v xl 
v x + 5v x2 
v x + 5v x3 



) 




(0 = Sv c 




(109) 
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where the shifts Sv v i, 5v v2 and 5v^ have been absorbed into the redefinition of the unde- 
termined parameters v^, v v and respectively. The new vacua is obtained by searching 
for the zeros of the F-terms, i.e. the first derivative of w v + 5w v with respect to the driving 
fields x° ; V^ ; 0°; P° an d A . By keeping only the terms linear in the shift 5v and neglecting 
the terms proportional to 5v/A, the minimization equations become 

fi[v<p(6v x2 - 5v x3 ) - v^v^} + f2V x dv c + oiu x ^t;a/A = 

fi[v v (6v xl - Sv x2 ) - v x Sv v2 ] + f2V x 8v{ + a 2 v x v (t> v A /A = 

h[v v (5v x a - Sv x x) - v x Sv v2 \ + f 2 v x 5v c + a 3 v x v^v A /A = 

2f 3 v x (5v xl + 5v x2 + 5v x3 ) + favySvc, + bxv x v^v A /A = 

2f 3 v x (5v xl + 5v x2 + 5v x3 ) + 2f 4 v (p Sv ip2 - fsv^Sv^ + b 2 v x v^,v A /A = (110) 

where the coefficients 01,2,3 an d 61,2 are linear combinations of the subleading coefficients 

at = fa + k 2 + k 3 - (fa + fa)v tp /v x - fav^/vA + fav^/ (v x va) + (-fa + fa)v v /v<p 

+(k w - fax)v v v v / O x tV>) + kmVyVt/fov/s) 
a 2 = fa - 2k 3 + (fa - fa)v v /v x - fav^l (v x v A ) - 2fav v /v (j> - (fa + fa^v^/ (v x v^) 

+k 12 v v vt/(v (t> VA) 

a 3 = fa - fa + fa + fav^/v A + (fa + fa)vr,/v^ + k 1 2V r ,v^/(v^v A ) 

6l = -Wi-W 2 + W 5 v^/VA + (WQ-W 7 )v v /V ( p 

b 2 = —Wi +W 2 + W 4 V ip /v x + W 5 V^/VA + (w e + W 7 )v r] /v<p + W&VyVqVz/ '(v^va) (HI) 

The solution to the linear equations Eq. flllOp is 

r _ % - a 2 v x v (t> VA (ax + a 2 + a 3 )/ 4 /5 v^va (bx + h)hh ~ 2bxf2f4, v$v A 

5v x2 
Sv x3 
5v V 2 



3/i Av v 


18(/i/ B - hh)h 


Av x 


12(/i/ B - 


f2h)h A 


a 2 - a ± Vy/Vcj/VA 


(ax + a 2 + a 3 )f 4 f 5 


V^V^Va 


(bx + b 2 )fxf 5 


- 261/2/4 V^Va 


3/i Av v 


18(/l/ 5 -/ 2 /4)/3 


Av x 


12(/i/ B - 


/ 2 / 4 )/ 3 A 


ax - a 3 v x VfpVA 


(ax + a 2 + a 3 )f 4 f 5 


V^Va 


(bx + b 2 )fxf 5 


- 261/2/4 V^Va 



3/i At V 18(/i/ 5 - / 2 / 4 )/ 3 A^ x 12(/,./-, ./,/,)/:; A 

(a x + a 2 + a 3 )/ 5 i^ua (61 - 6 2 )/ 2 w x ^a 



3(/i/ 5 -/2/ 4 ) A 2(^/5-/2/4) At- 



_ (ai + a 2 + a 3 )/4^A (61 - 6 2 )/i v x v^va , n2 ^ 
3(^/5-/2/4) A 2(^/5-/2/4) A^ 1 j 

In the same way, we obtain the minimization equations for the shifts 5f ^1,2,3 and fo^i 

(-2gxi^ + g 2 v v )5v (j>3 + g 2 v^5v vl + cxv x v\j 'A 2 = 
(45-1^ + g^Sv^ + c 2 v 2 x v\jA 2 = 
(-20iU0 + g^Sv^x + c 3 f x f J/A 2 = 

2g 3 v cj) 5v <t>3 + 2g i v v Sv vl + c 4 w 2 ?;J/A 2 = (113) 
where the parameters 01,2,3,4 are given by 
ci = 6s 2 + 3(s 7 + ss)^/ + 2si 4 w 2 / u 2 + (sxe + 8x7)1^/ (v^v^) + (2s 19 - 4s 20 )tv/ u x 
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-2(s 23 + S24)1VV|,/(W) + 2s28V v v v /(v x v^ 



c 2 = 6si + 3s 6 f v / Vtf, + 4si 3 ^/^ + 2s 15 vlv v / (v x v^) + (2si 9 - 4s 20 )%/ ^ 

.2 //„, „,2> 



-2(s 23 - su^vj {v x v^) + 2s 28 v tp v^/ (v x v^) 
c 3 = 6s 2 + 3(s 7 - ss)v v /v^ + 2si 4 ^/^ + (sie - sv^v^/^v^) + (2s i9 + 8s 20 )^/f x 

c 4 = 3r 2 + 3r 6 uJ/t;J + r 8 uj/uj + ri ujuj/ (uju J) + 4r lx v v /v x + 2r 12 v ip v r) / {v x v^) (114) 
Obviously Eg. (I113j) admits the solutions 



5^2 



c 3 ffig 
4<?i A 2 

c 2 vjjity 



2^ A 2 

16^4-2^3 A 2 
_ c x g 2 g 3 + 2c^g x g 2 v\v< i> 
Sgfgi-gjgs A 2 



5U03 



Ki = -TZ 7T" ( 115 ) 



From the above equations Eg. (11 151) . we can clearly see that all the shifts Sv^i/v^, Sv^/v^, 
dvffi/Vfj, and Sv v i/v v are of order \\. Finally the equations for the corrections 5i>ai,2,3 are 

{Ah x v A + h 2 v $ )Sv A1 + d lV 2 v 2 JA 2 = 

2„,2 /a2 



-2hiV A + h 2 v^)5v A3 + d 2 v x v A /A =0 



x 

,2„,2 



2M A + h 2 v^)6v A2 + d 3 v 2 v 2 A /A 2 = (116) 



where the coefficients di 2 3 are 



cZi = 6*i + 3t 6 v^/v A + (2tio + 8ti 2 )v ip /v x + At^v^vJ (v x v A ) + 2t i7 v^vj/ (v x v A ) 
d 2 = 6t 2 + 3t 7 v^/v A + (2ti - At 12 )v <fi /v x - 2*i 5 t>^/ (u x u A ) + 2i 17 ^u|/ (%u A ) 



J/^ + 2* 20 ^/(^ A ) 

2 /„.2 , + „.2„. /^„.2„ 

2 //„. „.2 



•2t 19 v v /v x + t 21 v^/(v x v A ) 



d 3 = 6* 2 + 3t 7 v^/v A + (2*io - 4*i 2 )^/f x - 2t 15 v v v^/ (v x v A ) + 2tx 7 v^/(v x v A ) 
i/vl + hiv^/ivl 

The solutions to the above equations Eg. (I116j) are given by 



+2t 19 v 2 Jv 2 + t 21 v 2 v<:/(v 2 x v A ) (117) 



di v 2 v A 
5VA1 = -2h^ 



'A2 



4/ii A 2 



It is obvious that Sv A i t2t3 /v A are of order Af, this is because the corrections to the vacuum 
alignment of A and £ arise at the next to next leading order. In short summary, the 
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modified vacuum configuration of the flavon fields can be parameterized by Eq. (11091) . 
the shifts Sv^i, 5v v2 and 5v% have been reabsorbed into the redefinition of v v , v v and 
respectively, which remain undetermined. The subleading corrections are suppressed by at 
least one power of 1/A with respect to the LO results, concretely Sv x i s 2,3/v x , Sv^/v^ and 
8v^/v x are of order A 2 , while Sv^x^/v^, 8v v ijv n and dv^i^/v^ are of order A^. These 
order of magnitudes can be clearly seen from Eqs. (II 121)1 15111 18j) . we note that the different 
suppressions of the shifts are due to the constraint of the flavor symmetry S A x Z% x Z A . 

Appendix C: GUT symmetry breaking 

In the following, we shall briefly discuss the GUT Higgs sector of the model in the 
present effective theory. Our Higgs sector is composed of H 5 , H A5 , H^, H-^ and if 24, the 
LO S4 x Z 3 x Z A invariant interactions between the different Higgs chiral superfields in 
the model are 

3 3 1 

w H = m 24 H 24 H 24 + X 24 H 24 H 24 H 24 + ^ fm—H^Of + ^ Xm—H^H^O^ 

i i 
2 1 1 3 1 

+ E cm^H 2A H 45 0f } + b' u —H M H 2A H^ + ]T b m —H M H 2i H,0 { P 

i i 

+ m 45 H m Hte + a H H M H Ab H 2A (119) 

where 



(7) = {X 3 ,XV,^} (120) 
Using the vacuum alignment shown in Eq. (ll02j) we can immediately obtain that 

wh = m 2A H 2A H 2A + \ 2 aM 2 aM 2A H 2A + fnH-^H 5 + X H H^H 5 H 24: + chH-^H^H^ 
+ bnH^H 24: H 5 + m&H^H 45 + anH^H^H^ (121) 



with 



f 9 f v\v x vav x v^ 

JH - Z JH1—^- + JH3 A2 

\ H - 2A m -^- + A^ 3 A3 

2 3 

b H = Qb H2 ^f + 2b H3 ^ (122) 



Since all the Higgs fields are neutral under the continuous U(1)r symmetry, the superpo- 
tential Eq. (1119p explicitly break U(1)r, while preserve the usual R-parity. Certainly we 
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can construct invariant operators comprising the driving fields, the Higgs fields and an 
arbitrary number of flavon fields, however, these operators don't contribute to the scalar 
potential due to the vanishing VEVs of the driving fields. Consequently, To completely 
understand the GUT symmetry breaking, maybe we should go beyond the effective theory 
framework and consider the ultraviolet completion Q. We note that the effective superpo- 
tential in Eq. (II 19[) could help us to qualitatively understand the GUT symmetry breaking, 
although this approach is not so satisfactory because of the U(1)r symmetry breaking. 
In the context of the ultraviolet completion of the effective model, the terms in Eq . (1 1 1 9 p 
could be generated from a U(1)r conserving superpotential in which the breaking is me- 
diated by additional fields which carry U(1)r charge. The ultraviolet completion of the 
model deserves considerable theoretical work (please see Ref. [5H] for an example of the 
ultraviolet completion of the A4 model), it is beyond the scope of the present work. 

The scalar potential of the model is determined by the SUSY F terms, D terms and 
soft terms contributions. We notice that the first two terms in Eq. (I122p is the interactions 
for H 2 4, they are exactly the same as those in the conventional GUT theory, this is because 
that the Higgs if 24 is neutral under the flavor symmetries, consequently the SU (5) GUT 
symmetry is broken into the standard model one as usual. Subsequently the VEVs of H 5 , 
f/45, H$ and H-^ break the standard model symmetry into the residual SU(3) C x U(l) em . 
Recalling that the parameter tan /3 could be small or large in the minimal supersymmetric 
standard model, this means that a hierarchy between the Higgs VEVs v u and Vd can 
be accommodated. In exactly the same way, the minor hierarchy between v§ and v^g in 
Eq. (1351) can be achieved by moderately fine-tuning the parameters in the superpotential 



7 The same is true for a large class of models with discrete flavor symmetry, where a continuous U(1)r 
symmetry is used to solve the vacuum alignment problem. 
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